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1. Recipe for the Universe 

All the popular models for structure formation are based on three key 
ingredients: (a) a model for the background universe (b) some mechanism 
for generating small perturbations in the early universe and (c) specification 
of the nature of the dark matter. 

The background universe is usually taken to be a Friedmann model 
with an expansion factor a{t). Such a model is completely specified if the 
composition of the energy density and the Hubble constant are specified. 
We will take Hq = WOh kms^^ Mpc^^ and express the energy density 
of the various species in terms of the critical density pc = {SHq /SnG) = 
1.88/i^ X 10~^^gcm~^, by writing pi = QiPc for the i*^ species. From various 
observations, we can impose the following constraints: (i) 0.011/i~^ ^ ^ 
0.016/1-2 (ii) ^ 0.8 (iii) Qi^^ ~ 0.007/i (iv) = A.85h-^ x lO^^ (v) 
Qiotgj = 0, ^ 0.3. Theoretical models strongly favour $7 = 1 and it is usual 
to invoke either a cosmological constant and/or nonbaryonic dark matter 
to achieve this[2]. We shall denote by 0,dm the total contribution due to 
all nonbaryonic energy densities. 

Models for structure formation also need to assume that small pertur- 
bations in the energy density existed at very early epochs. These pertur- 
bations can then grow via gravitational instability leading to the struc- 
tures we see today. In most of the models these perturbations are gen- 
erated by processes which are supposed to have taken place in the very 
early universe (say, at z ^ 10^^). Infiationary models - which are probably 
the most successful ones in this regard - can produce density perturba- 
tions with an initial power spectrum Pin{k) — Ak. Since each logarith- 
mic interval in k space will contribute to the energy density an amount 
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Alik) = da^/d{lnk) = {k^P{k) /2ti'^) we find that oc k^ for P oc fc. 
The contribution to gravitational potential from the same range will be 

= A^(9i?o/4/c''a^) which is independent of /c if A^ oc /c^. Such a "scale- 
invariant" spectrum is produced in some other seeded models as well. All 
these models need to be fine-tuned to keep the amplitude A of the fluctu- 
ations small upto, say, z ^ 10^. 

Given a Friedmann model with small inhomogeneitics described by a 
power spectrum P{k, zin) at a high redshift z = z\^, we can predict unam- 
bigiously the power spectrum P{k, zd) at the epoch of decoupling zd ~ 10^. 
This is because the perturbations at all relevant scales are small ak, z zd 
and we can use linear perturbation theory during this epoch. The shape 
of the spectrum aX z = zd will not be a pure power law since the gravita- 
tional amplifiction is wavelength-dependent. In general, the power at small 
scales is suppressed (relative to that at large scales) due to various physical 
processes and the exact shape of the spectrum sX z = z-q depends on the 
kind of dark matter present in the universe. In a universe dominated by 
"hot dark matter" particles of mass m ~ 30ey, the power per logarithmic 
interval in /c — space, A(A;) = {k^ P{k) /2tt'^Y^'^ , is peaked at A; = /cmax = 0.11 
Mpc~^(m/30ey) and falls exponentially for k > fcmax- Hence, in these mod- 
els, the scale k = /Cmax will go nonlinear first and smaller structures have to 
form by fragmentation. If the universe is dominated by "cold dark matter" 
particles with mass m ^ 35GeV, then A(A;) is a gently increasing func- 
tion of k for small k. If we set P{k) oc fc" locally, the index n changes 
from 1 at ^ 200/i~^Mpc to at ~ 10/i~^Mpc and to about 
(—2) at k~^ ~ lh~^ Mpc. In such models small scales will go nonlinear 
first and the structure will develop heirarchically[8]. 

The situation is more complicated if two different kinds of dark matter 
are present or if the cosmological constant is nonzero. The presence of 
the cosmological constant adds to the power at large scales but suppresses 
the growth of perturbations at small scales. Similar effect takes place if 
a small fraction of the dark matter is hot and the bulk of it is cold (eg. 
^^HDM — 0.2, ^cDM — 0.8). In both the cases there will be more power 
at large scales and less power at small scales, compared to standard CDM 
model. The spectrum A(fe) is still a gently increasing function of k and 
small scales go nonlinear first. 

The fact, that one can compute the power spectrum at z ^ Z£, analyt- 
ically, allows one to predict large scale anisotropics in CMBR unambigu- 
ously in any given model. Comparing this prediction with the anisotropy 
observed by COBE one can fix the amplitude A of the power spectrum. 
For a wide class[7] of the models, A(A:) = lQ-^{kLf with L ~ (24 ± A)h-'^ 
Mpc for A;-^ ^ SO/i"^ Mpc. For CDM like models the function A{k) flat- 
tens out at larger k and is about unity around k~^ ~ 8/i~^ Mpc. In pure 
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HDM models, A(/c) has a maximum value of ~ 0A2h-^im/30eVf at 
km — 0.11Mpc~^(m/30ey) and decreases exponentially at A; ^ km- 

The evolution of the power spectrum after decoupling (for z < zd) is 
more difficult to work out theoretically. In general, the power spectrum 
grows in amplitude (preserving the shape), as long as the perturbations 
are small[8]. In this case, we can write A{k,z) = [f{z)/f{zD)] A{k,ZD) for 
z < ZD- For example, in CDM models with = l,f{z) = (1 + z)~^; thus 
A (A;) grows by a factor 10^ at all scales between the epoch of decoupling 
{zd — 10^) and the present epoch (z = 0), if we assume that linear theory is 
valid at all scales. The resulting Ao(A;), obtained by linear extrapolation, is 
often used to specify the properties of the models. This spectrum correctly 
describes the power at large scales (say, k~^ ^ 30/i~^ Mpc) where Aq ^ 0.1. 
The "density contrast" cr(i?) measures the rms fluctuations in mass within 
a randomly placed sphere of radius R; upto factors of order unity, o"(i?) ~ 
A(A; ~ R~^) in heirarchical models. For most of the, COBE normalised, 
CDM-like models a{R) k, 1 around R 8h~^Mpc. Clearly linear theory 
cannot be trusted at smaller scales. 

There arc two major difficulties in understanding the physics at these 
small scales. Firstly, the true power Atrue(^) of dark matter will be larger 
than Ao(A:) due to nonlinear effects which are difficult to model analytically. 
Since dark matter particles interacts only through gravity, it is, of course, 
possible to study the formation of dark matter structures by numerical 
simulations. But to gain insight into the dynamics, it will be helpful to 
have simple analytic models explaining the N-body results. 

Secondly, it is important to understand gas dynamical processes before 
one can compare theory and observations at small scales. Since baryons 
can dissipate energy and sink to the minima of the dark matter potential 
wells, the statistical properties of visible galaxies and dark matter halos 
could be quite different. The situation is further complicated by the fact 
that in hierarchical models, considerable amount of merging takes place at 
small scales. It is usual to quantify our ignorance at these scales by a 'bias' 
(acronym for 'Basic Ignorance of Astrophysical Scenarios') factor b and 
write S.gai{r) = b"^ S,mass{^) - Such a parametrisation is useful only if b is ide- 
pendent of scale and morphology of galaxies. This seems to be somewhat 
unlikely. Since small scale observations are based on galactic properties, 
while theoretical calculations usually deal with underlying mass distribu- 
tion, any scale (or morphology) dependence of b could play havoc with 
predictive power of the theory. 

Recently some amount of progress has been achieved as regards the 
first aspect, viz, understanding nonlinear clustering of dark matter [6]. This 
approach is based on the relationship between the mean correlation function 
^{x,a) and the mean relative pair velocity v{x,a). These quantities are 
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related by an exact equation. 

dF dF 

where F = In + ^)] , A = hia,X = hix and h = —{v/ax). The char- 
acteristics of this equation shows that, as the evolution proceeds, power 
from a large scale I is transferred to smaller scales upto x = 1{1 + ■ 
By analysing the behaviour of h, it is possible to express ^(x,a) in terms 
of the mean correlation function in the linear theory^ ^i(/,a). It turns out 
that: i{a,x) = Q[h{a,l)T with = x'^{l + i) where Q = n = 1 for 
ih < 1-2; Q = 0.7, n = 3 for 1.2 ^ < 6.5 and Q = 11.7, n = 1.5 for 

> 6.5. This relation shows that ^ is steeper than ^l. 

Unfortunately, no such simple pattern exists in the dynamics of baryons 
coupled to dark matter. The gas dynamical processes introduce several 
characteristic scales into the problem and the evolution becomes quite com- 
plicated. The only reliable way of probing these systems seems to be through 
massive hydro simulations which are still at infancy. 

It is clear from the above discussions that our theoretical understand- 
ing is best at large scales {k~^ ^ 30h~^ Mpc) where linear theory is valid, 
Ao(A;) is well determined and baryonic astrophysical processes are not im- 
portant. At the intermediate scales (3/i~^ Mpc^ ^ 30h~^ Mpc), it is 
not too difficult to understand the dark matter dynamics by some approx- 
imation but the baryonic physics begins to be nontrivial. At still smaller 
scales, {k~^ ^ 3h~^ Mpc) there is considerable uncertainty in our theoreti- 
cal predictions. We shall now turn to the observational probes of the power 
spectrum at different scales. 



2. Probing the power spectrum 

One of the direct ways of constraining the models is to estimate the density 
contrast f7obs(^) from observations at different scales and compare it with 
the theoretically predicted values. FortTinatcly, we now have observational 
probes covering four decades of scales from 10^^ Mpc to 10"^ Mpc. We shall 
discuss the probes of different scales in the decreasing order. 



2.1. NEAR HORIZON SCALES: (300 - 3000) H-^MPC 

These scales are so large that the best way to probe them is by studying the 
MBR anisotropy at angular scales which correspond to these linear scales. 
Since a scale L subtends an angle 9{L) = 1°{L/I00h~^ Mpc) at 2; ~ zd, 
the (AT/T) observations at (3° - 30°) probe these scales. The COBE- 
DMR observations[ll] of (AT/T)rms and {AT/T)q allow one to obtain the 
following conclusions: (i) a{10^h~^ Mpc)~ 5x 10~^ (ii) The power spectrum 
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at large scales is consistent with Pm{k) Ak and, if we take = 1,, then 
^1/4 ^ (24 ± 4)h-^ Mpc (iii) In this range, a{R) ^ (24 ± Ah'^ Mpc/Rf. 

2.2. VERY LARGE SCALES : (80 - 300)H-^ MPC 

2.2.1. CMBR probes: 

These scales span (0.8° — 3°) in the sky at z ~ zd- Several ground based and 
balloon-borne experiments to detect anisotropy in MBR probe this scale. 
For example, the UCSB South Pole experiment has reported [10] a prelim- 
inary 'detection' of (AT/T) ~ 10"^ at 1.5° scale, and a 95% confidence 
level bound of (AT/T) < 5 x 10^^. This translates into the constraint of 
(7(102/i-i Mpc) ^ 5 X 10-2. 

The angular anisotropy of CMBR is dominated by the gravitational po- 
tential wells of dark matter at large scales. However, at ^ ~ 1° baryonic 
process affect the pattern of anisotropy significantly. The precise deter- 
mination of degree scale anisotropy can, therefore, help in distinguishing 
between different models [14]. 

2.2.2. Galaxy surveys: 

Some galax;y surveys, notably CfA2 survey and pencil-beam surveys probe 

scales which arc about 10^h~^ Mpc in dcpth[l]. Unfortunately, the statistics 
at these large scales is not good enough for one to obtain cr{R) directly from 
these surveys. 

2.3. LARGE SCALES : (40 - 80) R-^ MPC 

2.3.1. CMBR probes: 

The scales correspond to Ombr — (24' — 48') and are probed by the ex- 
periments looking for small angle anisotropics in MBR. The claimed de- 
tection[3] by MIT-MASM of (AT/T) ^ (0.5 - 1.9) x 10"^ at 61 ~ 28', if 
confirmed, will give a bound of (7(50/i-^ Mpc ) 0.3. 

2.3.2. Galaxy Surveys: 

Several galaxy surveys, in particular the IRAS-QDOT and APM surveys, 
give valuable information about this range [9]. The angular correlation of 
galaxies, measured by APM survey is u}{d) ~ (1 — 5) x 10"^ at ^ ~ 14°. This 
corresponds to a{^{)h~^ Mpc) = 0.2. What is more important, these surveys 
can provide valuable information about the shape of the power spectrum in 
this range if we assume that galaxies faithfully trace the underlying mass 
distribution. 
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2.3.3. Large scale velocity field: 

Using distance indicators which are independent of Hubble constant, it 

is possible to determine the peculiar velocity field v{R) of galaxies upto 

about 80h-^ Mpc or so. The motion of these galaxies can be used to map 

the underlying gravitational potential at these scales. Careful analysis of 

observational data shows[5] that v{AOh~^ Mpc) ~ (388 ± 67) kms~-^ and 

v{60h~^ Mpc) ~ (327ib82) kms^-*^. From these values it is possible to deduce 

that a{50h~^ Mpc ) ~ 0.2. These observations also allow us to determine 

the value of the parameter (fi°'^/6iRAs) where 6iras is the bias factor with 

respect to IRAS galaxies. One finds that (r^°-^/6iRAs) = l-28lo:59 which 

implies that if O = 1, then 6iras = 0-781q 29 ^IRAS = 1 then Q = 

1 ^1+1-74 

^■'^^-0.97- 

2.3.4. Clusters and voids: 

The cluster-cluster corelation function and the spectrum of voids in the 
universe can, in principle, tell us something about these scales. Unfortu- 
nately, the observational uncertainties are so large that one cannot yet 
make quantitative predictions. 



2.4. INTERMEDIATE SCALES : (8 - 40)^?-^ MPC 

2.4.1. Galaxy Surveys: 

The galaxy - galaxy correlation function ^gg = [r /5Ah~^ Mpc]~^'^ is fairly 
well determined at these scales. Direct observations suggest that crgai(8^~^ 
Mpc) ~ 1 but the (Tdm and CTgai at these scales can be quite different because 
of possible biasing. 

2.4.2. Cluster Surveys: 

There have been several attempts to determine the correlation function of 
clusters of different classes. It is generally believed that ^cc — {r/L)~^'^ 
with L ~ 25h^^ Mpc. The index n = 1.8 is fairly well determined though 
the scale L is not; in fact, L seems to depend on the richness class of 
the cluster. The quantity (Ccc/Cgg)"^''^ can be thought of as measure of the 
relative bias between cluster and galaxy scales. Observations suggest [4] that 
this quantity depends on the cluster class and varies in the range (2 — 8). 
The observational uncertainties are still quite large for this quantity to be of 
real use; but if the observations improve we will have valuable information 
from ^cc- 

2.4.3. Abundance of rich clusters: 

The scale R = Sh'^ Mpc contain a mass of 1.2 x lO^^Qh^QMQ. When 
this scale becomes nonlinear, it will reach an overdensity of about S ~ 
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178, or - equivalently - it will contract to a radius of Rf ~ {8h ^ Mpc) 
/(178)i/3 ~ 1.5/1-1 Mpc. A mass of IO^^Mq in a radius of 1.5 Mpc is a good 
representation of Abell clusters we see in the universe. This implies that 
the observed abundance of Abell clusters can be directly rclaied to a{8h^^ 
Mpc). Several people have attempted to do this [13]; the final results vary 
depending on the modelling of Abell clusters, and give a{8h''^ Mpc) ~ 
(0.5-0.7). Since crgax{8h~'^ Mpc)~ 1, this shows that b ~ (1.23-2) at Sh'^ 
Mpc. 

It is possible to give this argument in a more general context[12]. Sup- 
pose that the contribution to critical density from collapsed structures with 
mass larger than M is U{M), at a given redshift z. Then one can show that 



n{M) = erfc 



V2ao{M) 



where 5c = 1.68 and erfc(x) is the complementary error function. The Abell 
clusters (at z = 0) contribute in the range CI 2± (0.001 — 0.02). Even with 
such a wide uncertainty, we get adus — (0-5 — 0.7). 

2.5. SMALL SCALES : (0.05 - 8)iJ-i MPC 

These scales correspond to structures with Msmooth ~ (3 x 10^ — 1.2 x 
10^5) nh^gMQ and we have considerable amount of observational data 
covering these scales. Unfortunately, it is not easy to make theoretical pre- 
dictions at these scales because of nonlinear, gas dynamical, effects. 

2.5.1. Epoch of galaxy formation: 

Observations indicate that galaxy-like structures have existed even at z ~ 3. 
This suggests that there must have been sufficient power at small scales to 
initiate galaxy formation at these high rcdshifts. Unfortunately, we do not 
have reliable estimate for the abundance of these objects at these redshifts 
and hence we cannot directly use it to constrain a{R). 

2.5.2. Abundance of quasars: 

The luminosity function of quasars is fairly well determined upto z ~ 4. If 
the astrophysical processes leading to quasar formation are known, then the 
luminosity function can be used to estimate the abundance of host objects 
at these redshifts. Though these processes are somewhat uncertain, most 
of the models for quasar formation suggest that we must have a{0.5h~^ 
Mpc) ^ 3. 
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2.5.3. Absorption systems: 

The universe at 1 ^ 2; ^ 5 is also probed by the absorption of quasar hght 
by intervening objects. These observations suggest that there exist signif- 
icant amounts of clumped material in the universe at these redshifts with 
neutral hydrogen column densities of A^hi — (10^^ — 10^^)cm~^. We can con- 
vert these numbers into abundances of dark matter halos by making some 
assumptions about this structure. We find that [12] in the redshift range 
of 2; ~ (1.7 — 3.5) damped Lyman alpha systems contribute a fractional 
density oi ^Ly ^ (0.06 - 0.23). This would require ^(lO^^Mo) ~ (3 - 4.5). 

2.5.4. Gunn-Peterson bound: 

While we do see absorption due to clumped neutral hydrogen, quasar spec- 
tra do not show any absorption due to smoothly distributed neutral hydro- 
gen. Since the universe became neutral at 2; ^ z/j ~ 10^, and since galaxy 
formation could not have made all the neutral hydrogen into clumps, we 
expect the IGM to have been ionised sometime during 5 ^ z ^ 10^. It 
is not clear what is the source for these ionising photons. Several possible 
scenarios (quasars, massive primordial stars, decaying particles etc.) have 
been suggested in the literature though none of these appears to be com- 
pletely satisfactory. In all these scenarios, it is necessary to form structures 
at z ^ 5 so that an ionising flux of about J = lO^^^ergs cm^^s^^ Hz^^ sr~^ 
can be generated at these epochs. Once again, it is difficult to convert this 
constraint into a firm bound on a though it seems that a{0.5h-^ Mpc) ^ 3 
will be necessary. 

3. Gravitational lensing and large scale structure 

In the above discussion we have not taken into consideration the constraints 
imposed by gravitational lensing effects on the structure formation models. 
This aspect will be discussed in detail in the other articles in this volume; 
here we shall contend ourselves with a brief mention of the possibilities. 

Gravitational lensing probes the gravitational potential directly and can 
provide valuable information at very different scales. At the largest scales 
{R ~ 10^ Mpc) lensing can be used to probe the geometry of the universe. 
For example, it is possible to put firm bounds on the energy contributed 
by cosmological constant from such considerations. 

At intermediate scales ( R ~ 50Mpc) lensing has the potential of pro- 
viding information about the power spectrum of fluctuations which are in 
the quasilinear phase. In principle the distortion of images can be inverted 
to obtain this information, though in practice this is extremely difficult. 

At smaller scales, the "weak lensing" - leading to arcs and arclets at 
cluster scales - is already providing a clue to the mapping of dark matter 
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distribution in clusters. On the other hand, direct optical and X-ray ob- 
servations provide us information about the distribution of visible matter 
in clusters. The combination of these techniques should give us valuable 
information as regards the dynamical processes which separated baryons 
from dark matter. 

At still smaller scales, galactic potentials have the capacity to produce 
multiple images of distant sources. The statistics of these multiple images 
depends crucially on the core radii of the galaxies, which in turn depends 
sensitively on the structure formation models. The absence of significant 
number of multiple images with large angular separations puts severe con- 
straints on models for structure formation. The analytic modelling of non- 
linear dark matter clustering described earlier could be used to strengthen 
these constraints still further. 

4. Scorecard for the models 

The simplest models one can construct will contain a single component of 
dark matter, either cold or hot. Such models are ruled out by the observa- 
tions. The HDM models, normalised to COBE result will have maximum 
power of ^ 0A2h-'^{m/30eVf at k = km = 0.11 Mpc"^(m/30ey). In 
such a case, structures could have started forming only around {1 + Zc) = 
(Am/1.68) = h^Q{m/30eV)'^ or at Zc = 0. We cannot explain a host of 
high-z phenomena with these models. The pure CDM models face a differ- 
ent difficulty. These models, normalised to COBE, predict as — I, which is 
too high compared to the bounds from cluster abundance. When nonlinear 
effects are taken into account, one obtains ,^gg oc r~^'^ for h = 0.5 which is 
too steep compared to the observed value of ^gg oc r~^'^. In other words, 
CDM models have wrong shape for ^(r) to account for the observations. 

The comparison of CDM spectrum with observations suggests that we 
need more power at large scales and less power at small scales. This is 
precisely what happens in models with both hot and cold dark matter or 
in models with nonzero cosmological constant. These models have been 
extensively stTidicd during the last few years, and they fare well as far 
as large and intermediate scale observations are concerned. However, they 
have considerably less power at small scales compared to CDM model. 
As a result, they do face difficulties[12] in explaining the existence of high 
redshift objects like quasars. For example, a model with 30% HDM and 70% 
CDM will have (T0.5 1.5; to explain the abudnance of quasars comfortably 
one needs o"o.5 c± 3.0. To explain the abundance of damped Lyman alpha 
systems one requires still larger vlaues of about cro.5 4 or so. Demanding 
that C7(1O^^M0) > 2 [which is equivalent to saying that lO^^M© objects 
must have collapsed at a redshift of Z12 = (2/1.68) — 1 ~ 0.2] will completely 
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rule out this model. Similar difficulties exist in models with cosmological 
constant. Notice that all models are normalised using COBE results at very 
large scales. Hence the severest constraints arc provided by observations at 
smallest scales, since the "lever-arm" is longest in that case. 

The comparison of models show that it is not easy to accommodate all 
the observtions even by invoking two components to the energy density. 
(These models also suffer from serious problems of fine-tuning). By and 
large, the half-life of such quick-fix models seem to be about 2-3 years. 
One is forced to conclude that to make significant progress it is probably 
necessary to perform a careful, unprejudiced analysis of: (a) large scale ob- 
servational results and possible sources of error and (b) small scale baryonic 
astrophysical processes. 
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